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I. INTRODUCTION 



The observation of the by the LEPS collaboration [1] and its subsequent confirmation 
have brought a lot of excitements in the field of hadronic physics [2] . On the other hand, there 
are increasing number of experiments reporting negative results. In particular, the latest 
experiments at JLAB[3] find no signal from the photoproduction process on a deuteron nor 
on a proton target, from which the G"*" was observed earlier by the SAPHIR collaboration 
with lower statistics. Although the present experimental results are quite confusing and 
frustrating[4] , one can not afford to give up further refined experimental search, because if a 
pentaquark is found, it will provide a major and unique testing ground for QCD dynamics 
at low energy. 

Another multiplet to search for as possible pentaquark states are those with one heavy 
antiquark. The HI collaboration at HERA has recently reported on the finding of an anti- 
charmed pentaquark 0c(3O99) from the D*p invariance mass spectrum[5]. Unfortunately 
other experiments could not confirm the finding[6-8]. While the experimental search for 
the heavy pentaquark is as confusing as that for the light, theoretically, the heavy and 
light pentaquarks stand on quite different grounds. Cohen showed that the original pre- 
diction for the mass of the based on the SU(3) Skyrme model[9] is not valid because 
collective quantization of the model for the anti-decuplet states is inconsistent in the large 
Nc limit [10]. In contrast, many theories consistently predicted a stable heavy pentaquark 
state. The pentaquark with one heavy anti-quark was first studied in Ref. [11, 12] in a 
quark model with color spin interaction. Then it has been studied in quark models with 
flavor spin interaction[13] and Skyrme models [14, 15], and with the recent experiments, at- 
tracted renewed interests [16-18], some of which were motivated by the diquark-diquark[19] 
and diquark-triquark[20] picture. Such states also appear naturally in a coupled channel 
approach[21], and in the combined large Nc and heavy quark limit of QCD[22]. If the heavy 
pentaquark state is stable against strong decay, as was predicted in the D meson bound 
solition models [14], it could only be observed from the weak decay of the virtual D me- 
son. From a constituent quark model picture based on the color spin interaction [23], one 
expects a strong diquark correlation, from which one could have a stable diquark-diquark- 
antiquark[19] or diquark triquark[12] structure. The question is whether such strong diquark 
structure will survive other non-perturbative QCD dynamics in a multiquark environment 
and produce a stable pentaquark state. Such questions are being intensively pursued in 
quark model approaches [24-26]. In particular, an important question at hand is whether 
the net attraction from the diquark correlations in the pentaquark configuration is stronger 
than that from the corresponding diquark and additional quark-antiquark correlation present 
when the pentaquark separates into a nucleon and a meson state. Since the correlation are 
inversely proportional to the constituent quark masses involved, the attraction is expected 
to be more effective for pentaquark state with heavy antiquark. Another non-perturbative 
approach that can be used to answer such question is the QCD sum rule method. 

There have been several QCD sum rule calculations for the light pentaquark states[27-33]. 
The application to the heavy pentaquarks was performed by two of us in a previous work[34], 
where we used a pentaquark current composed of two diquarks and an antiquark, and found 
the sum rule to be consistent with a stable positive parity pentaquark state. The similar 
approach has been applied to the sum rules for Z}s(2317) [35]. In this work, we extend the 
previous QCD sum rule calculation to investigate the anti-charmed pentaquark state with 
and without strangeness using two different currents for each case. We find a convergent 
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Operator Product Expansion (OPE) only for the non-strange heavy pentaquark sum rule 
obtained with an interpolating field composed of two diquarks and one anti-charm quark, 
that has been previously used by us [34]. The stability of non-strange heavy pentaquark 
is consistent with the result based on the quark model with flavor spin interaction [25]. We 
then refine the convergent sum rule by explicitly including the DN two-particle irreducible 
contribution. The importance of subtracting out such two-particle irreducible contribution 
has been emphasized in Ref. [36-38] for the light pentaquark state. In fact, estimating the 
contribution from the lowest two-particle irreducible contribution is equally important in 
lattice gauge theory calculations [39, 40] to isolate the signal for the pentaquark state from 
the low-lying continuum state. We find that for the non-strange heavy pentaquark sum rule, 
including the DN continuum contribution tends to shift the position of the pentaquark state 
downwards. Given the negative experimental signatures of the charmed pentaquark states 
above the threshold, the present result suggests that the anti-charmed pentaquark states 
might be bound as was predicted in D meson bound soliton models. 

This paper is organized as follows. In Section 11, we introduce the interpolating field 
for the 0c and discuss the dispersion relations that we will be using. Section III gives the 
phenomenological side and Section IV gives the OPE side. The QCD sum rules for ©c and 
their analysis are given in Section V. 

II. QCD SUM RULES 

A. Interpolating field for ©c 

Let us introduce the following two interpolating field for ©c, 

©cl = e'''"'{u^C-ff,Ub)j5-f^dc{cdij5dd) , 

©,2 = e^''{e'^^^ulCj,df){e'^'^u^Cdn)Ccl . (1) 

Here the Roman indices a,b, . . . are color indices, C denotes charge conjugation, T transpose. 
Note that ©d is composed of a nucleon current (proton) and a pseudo scalar current (D), 
while ©c2 is composed of diquark-diquark-antiquark and has been investigated in a previous 
work [34] . 

For the charmed pentaquark with strangeness, we consider the following two possible 
currents, 

= e'''\e'''fu^Csf){e'^''u^Cd,)Ccl , 
©CS2 = e''''{e'''fulC^,Sf){e'^''u^Cd,)Ccl . (2) 

Here, instead of choosing Qcsi a-s a direct product of a nucleon and a Dg or a hyperon and a D 
meson currents as in ©d, we choose it to well represent a state having two diquark structure 
with the same scalar quantum number but with different fiavor. Such configuration allows 
all the five constituent quarks to be in the s-wave states, which will have the lowest orbital 
energy and consequently could be the dominant ground state configuration [18]. Moreover, 
as we will see, ©d couples dominantly to the nucleon and D meson state, suggesting that 
currents composed of a direct product of a nucleon and a meson currents are not suitable 
for investigating the properties of the pentaquark state. 
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Under parity transformation q'{x') — Joq{x), the ©c currents transform as, 

e>;.i = -ToOci , ©;,2 = 7oec.2. (3) 

B. Dispersion relation 

The first type of QCD sum rules for the heavy pentaquarks that we will be using are 
constructed from the following time ordered correlation function, 

UT(q)^i j d'xe"^-%Q\T[Q,{x),QMm = M^^) + i^.i^^) , (4) 

where Oc can be any of the currents in Eq. (1) or in Eq. (2), and 11^, IIi are called the 
chiral even and chiral odd parts respectively. As can be seen in Eq. (3), the currents are 
not eigenstates of the parity transformation and can couple to both positive and negative 
parity states. The spectral densities calculated from the OPE of Eq.(4) are matched to that 
obtained from the phenomenological assumption in the Borel-weighted dispersion integral, 

fSo 1 

/ dg^e-^^/^^V(g^)-im[nf ^"^(g^) - nr(g^)] = , = 1, g) , (5) 

where is the Borel mass. Here, higher resonance contributions are subtracted according 
to the QCD duality assumption, which introduces the continuum threshold Sq. We have 
also introduced an additional weight function W{q'^) for later use. 

In this work, we will also work with the "old-fashioned" correlation function, which is 
defined as [29] 

UT{q)^i J dW«-^(0|^(x°)ee(x)ee(0)|0) . (6) 

This type of correlation function has been used in projecting out positive and negative parity 
nucleon states [41]. We then divide the imaginary part into the following two parts, which 
are defined only for > 0, 

-lmU{qo) = A(go)7° + i^(?o) ■ (7) 

TT 

One should note that these can be identified with the imaginary part calculated from Eq.(4), 

MQo) = -linUq{qo)qo 

TT 

B{qo) = -lmni(go), (8) 

TT 

for go > 0. 

Now, depending on the parity of the current ©c in Eq.(3), one can extract the positive 
or negative-parity physical state only by either adding or subtracting A and B. That is, the 
spectral density for the positive and negative parity physical states will be as follows, 

^±(„ \ _ j MQo) T B{qo) For Gd, 0^,1 
P ^^0^ - 1 A{qo) ± B{qo) For 0^2, ©cs2 " 
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The sum rules are then obtained by again matching the corresponding spectral density from 
the OPE and phenomenological side, 

poo 

/ dqoe-^o/M^[p%M - Pope(?o)] = . (10) 

III. PHENOMENOLOGICAL SIDE 

A. 6ci,6csi 

For ©ci current, the interpolating field couples to a positive parity state as, 

(0|e,i(a;)|e,(p) : P = +) = A+,ei 75f/e(p)e-'^-" , (11) 

and to a negative parity state as, 

(0|eei(x)|ee(p) : P = -) = A_,,iC/G(p)e-^^-^ . (12) 

Here, X±^ci denotes the couphng strength between the interpolating field and the physical 
state with the specified parity. Similar relations will hold for Qcsi- Using these, we obtain 
the phenomenological side of Eq. (4) separated into chiral even (11^) and odd (Hi) parts, 
which are defined to be the parts proportional to ^ and 1 respectively. 

As was first pointed out in Ref.[36], the correlation function can also couple to the DN 
continuum state, whose threshold could be lower than the expected ©c mass. Its phenomeno- 
logical contribution can be estimated by using, 

{0\&,i\DN{p)) = tXDN,ciUN{p) . (13) 
Combining these two contributions, we find 

nT,cl W) = -I^±,cl| -o - ^\^DN,cl\ dp-. + (14) 



q — rriQ J — [p — qf — mjj 

where the minus (plus) sign in front of me is for positive (negative) parity. The dots denote 
higher resonance contributions that should be parameterized according to QCD duality. It 
should be noted however that higher resonances with different parities contribute differently 
to the chiral-even and chiral odd parts [42] . Thus, 11^'^'^'^ and H^^'^^ constitute separate sum 
rules. For ©csi, the D meson should be replaced by the Dg meson. 

The corresponding spectral density for the pole and DN contributions are given respec- 
tively by 

-ImnP°J^^(g) = i\X±,ci\'^S{q^ -ml)T'me\\±,ci\^S{q'^ -ml), 

TT ' 



+1 Wil'J^ - V(mi, + ml) + {ml - m\)\ (15) 



We notice that the chiral-odd part has opposite sign depending on the parity while the chiral 
even part has positive-definite coefficient. 
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B. @c2,&cs2 

As can be seen in Eq.(3), Qc2 transforms differently compared to 0ci under parity. Thus, 
tfie couplings to the interpolating field are 

(O|ee2(:r)|0e(p) : P = +) = A+,,2f/e(p)e-*"" , 

{0\e,2{x)\ecip) : P = -) = X-,c2 75t/G(p)e-^^-^ . (16) 
Similarly, the couphng to the DN continuum state changes as follows, 

(0|ee2|i^iV(p)) = A^iv,c275t/iv(p). (17) 

Combining these changes, we find. 

Consequently, the spectral densities are, 

-lmU^''^^{q) = i\X±,c2\'^S{q^ - m%) ± me\\±^c2?8{q^ - m|), 

TT ' 

llmn?J,(g) = i\\nN,c2? '^' ^^^^ 'f"" <l' " + mf,) + (m^, - mlY 



-lA 



DiV,c2 



32 TT^ g4 

- 2g^(m^ + m^) + (m^ - mlf. (19) 



C. Phenomenological side 

The final form for the phenomenological side to be used in Eq.(lO) can be obtained from 
combining Eq.(15) or Eq.(19) according to Eq.(9), both of which are given in the following 
form, 

Pphen(5o) = |A±p5(go - me) + - qQ)pDN{(lo) + ^(90 - \fs^)ptonMo) , (20) 

where the usual duality assumption has been used to represent the higher resonance contri- 
bution above the continuum threshold a/sq; i.e., pfontiQo) = Pope(?o)- The spectral density 
for the two-particle irreducible part is given by 



PDN,ciilo) = ^^32^2^ \/iQo-mDy -m%^J{qo + 



^{qo±m^r-ml 



Pdn,c2{Qo) = *^3^f * ^ (go - ruD)^ - mjf^ (go + mo)^ - m% 
^{qoTm^r-ml 

We substitute the above into the Borel transformed dispersion relation in Eq.(lO). 
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u u 

u c z ^ d i > M ^ c r :r > 




FIG. 1: Schematic OPE diagrams for the current 0ci in Eq.(l). Each label corresponds to that 
in Eq.(24). The solid lines denote quark (or anti-charm quark) propagators and the dashed lines 
are for gluon. The crosses denote the quark condensate, and the crosses with circle represent the 
mixed quark gluon condensates, (c) represents diagrams proportional to gluon condensate with 
gluons lines attached to the light quarks only, (d) represents those where the gluons are attached 
to the heavy quarks only, while (e) represents those where one gluon is attached to the heavy quark 
and the other to a light quark in all possible ways, (f) and (g) represent all diagrams that contain 
the quark-gluon condensate. 

IV. OPE SIDE 
A. Gel 

Here, we present the result for ©d- To keep the charm quark mass finite, we use the 
momentum-space expression for the charm quark propagators. For the light quark part of the 
correlation function, we calculate in the coordinate-space, which is then Fourier-transformed 
to the momentum space in D-dimcnsion. The resulting light-quark part is combined with 
the charm-quark part before it is dimensionally regularized at D — 4. 

Our OPE is given by 

+n(-^) + n(») + n^*) + n« + n(^) , (23) 
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where the superscript indicates each diagram in Fig. 1. The imaginary part of each diagram 
is calculated as 



TT 



-Imn(")(g2) 



-Imn(^)(g2) 

TT 



TT 



-Imn(^)(g2) 



-Imn(^)(g2) 

TT 



TT 



TT 



-Imn('*)(g2) 



TT 



-Imn(*)(g2) 
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5!6!2i3 7r8 



1 



(l-u) 



-36u{l - u)[-L{u)f 



+ mq^^l-ufi-Liu)]"] +mJ2[-L{u)f} , 



4!4!297r8 



(l-u) 



fl6mcu[—L{u)Y 



du- 



-2>u{l-u)[-L{u)f 



3 • 4!2i%6 Jo (1 - uf 



+ ^me[-LH]3 



+4gV(l-«)2[-L(M)] 



+ mc 



11' 



4!5!3 • 2i%6 



du 



u 



{l-u) 



i (^96k(1 --u) + 5(1 -'u))[-L(m)]= 



-[-LH]=^ + 6gVl-t.)[-LH]2}, 



(g^g ■ Gq) 
3!4!2i07r6 



du 



{\-uf 



irric \^2u{l - u)[-L{u)Y 



-mu^[-L{u)f - 12(1 - u)[-L{u)f + 72g2-u(i - uf[-L{u)Y 
-'^[-L{u)f + ?,q^u\l-u)[-L{u)fY 



{qqf 



du- 



i I2u{l - u)[-L{u)f 



9 • 2%4 (1 _ uf 

-16gV(l - uf[-L{u)] + 3(1 - 'u)[-L(M)]2j + 21mc[-L{u)f 

^^dM{-^meM+ [-L(m)] -2g\(l-M)} , 



9 • 2%2 
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ilmn(^)(g2) = (-^ + 22me) 5{q^ - ml). (24) 

TT Zib 

Here the upper limit of the integrations is given by A = 1 — m1/(f and L{u) — (fu{l — 
u) — mlu. Our OPE calculation has been performed up to dimension 12 here. Up to 
dimension 5, we include all the gluonic contributions represented by the gluon condensate 
and the quark- gluon mixed condensate. Beyond the dimension 5, we have included only 
tree-graph contributions which are expected to be important among higher dimensional 
operators. Other diagrams containing gluon components are expected to be suppressed by 
the small QCD coupling. Therefore, the higher order tree-graphs, which are the higher order 
quark condensates, will be able to give us an estimate on how big the typical higher order 
corrections should be beyond dimension 5. The integrations can be done analytically but 
we skip the messy analytic expressions. For the charm-quark propagators with two gluons 
attached, we use the momentum-space expressions given in Ref. [43]. The Wilson coefficients 
for light-quark condensates come from (gg)", where n = 2, 3, 4. This is in contrast with the 
OPE for 0c2, where the Wilson coefficient are non-zero only for n = 4. 

The first important question to ask in the OPE is whether it is sensibly converging as an 
asymptotic expansion. For that, we choose to plot the Borel transformed OPE appearing in 
Eq.(lO) after subtracting out the continuum contribution, 

POO 

n(^) {M')= dqo e-^o/M^ [pji^) (go) - P^!^^ (^o)] = 0. (25) 

^0 

Here j = a, 6, c. denotes each contribution in the OPE in Eq.(24) after adding according to 
the rules in Eq.(9). 

We use the following QCD parameters in our sum rules [29, 44], 

rus = 0.12 GeV , = 1.26 GeV , 

(^^G^^ = (0.33 GeV)^ , {G^) = 0.045 GeV"^^ , 

(qq) = -(0.23 GeV)^ , (ss) = 0.8{qq) , 
{qga ■ Gq) = (0.8 GeV^) x {qq) , {sga ■ Gs) = (0.8 GeV^) x {ss). (26) 

Fig. (2) represents the OPE as defined in Eq.(25) with the imaginary part in Eq.(24). One 
notes that for the negative parity case, the perturbative contribution is only a small fraction 
of the OPE, and hence do not converge. For the positive parity case, the power corrections 
alternate in signs, and the gluon condensate, which represents the light diquark correlation, 
is only a small correction to the power correction. Hence, such structure, would hardly 
couple to a pentaquark state, and it is meaningless to perform a detailed QGD sum rule 
analysis. We present the result with the continuum threshold 5*0 = (3.3 GeV)^. This value is 
chosen in the range \/So = 3.2 — 3.6 GeV, which has been used to analyze the anticharmed- 
pentaquark sum rule in Ref [34] . However changing 5*0 does not change the relative strength 
of each contribution, and hence the conclusion of this section. We will therefore, analyze 
the subsequent OPE with the same threshold. 
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FIG. 2: OPE as defined in Eq.(25) for the current 9^ and So = (3.3 GeV)^ The left (right) figure 
is for positive (negative) parity case. The sohd Une (a) represents the perturbative contribution. 
The hne specified as OPE represents the sum of the power corrections only, (c) represents the 
gluon condensates. Other labels represent contribution from each term in Eq.(24). Here we plot 
only a few selected terms in the OPE. 
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The OPE for Qc2 are given in Ref.[34]. Here, we rewrite the result for completeness, 

r-A 1 



-Imn('*)(g2) 



IT 



-Imn(*)(g2) 



-Imn(^)(g2) 

TT 

TT 



du- 



5 ■ 5! 2i27r8 Jq (1 - u)5 

'3! 3! 2i07r6 (^^') / '^'^ {1 - uf 
'^{qQ^H + 5{q^ -m^ 



l-u)+m,} [-L{u)f 



1 — n) + iTLc} [—L(u) 



cl -I 



1 



du- 



u 



5! 3! 3-2i07r6 \ vr ' / Jq (1 - m) 
X { 3ml^{l -u)+ mdl - u){3 - 5u)q'^ + 2uml \ [-^(m)^^ 
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FIG. 3: Similar figure as Fig. 2 for the current Qc2- Here each label represents contribution from 
each term in Eq.(27). The gluon condensates (b) are the dominant power correction in the positive 
parity channel (left figure). 



1 



TT 



5! 4! 2V 



du- 



u 



+6mcg {2u — 1) — 2m 



l-u] 
3m + 1 



u 



5u 

Y ~ 



1 fl 



u 



3m 



+ 7 



(27) 



The diagrams corresponding to every term above, denoted by the superscripts (a) — (e), can 
be found in Ref. [34]. 

Fig. (3) represents the OPE as defined in Eq.(25) with the imaginary part in Eq.(27). 
As can be seen from the left figure, the OPE without the perturbative contribution is 
dominated by the gluon condensate coming from the light diquarks. This suggests that the 
diquark correlation is the dominant interaction among the quarks and heavy antiquark in the 
positive parity channel. Moreover, the perturbative contribution is larger than sum of the 
power corrections denoted as "OPE" in the figure. Therefore, the pentaquark could couple 
strongly to this current and a detailed QCD sum rule analysis is sensible. The situation 
changes for the negative channel, where the power corrections have alternating signs, and 
hence becomes less reliable. 
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FIG. 4: Schematic OPE diagrams for the currents Qcsi in Eq.(28) and @cs2 in Eq.(29). Each label 
corresponds to that in Eq.(28) or Eq.(29). All the other notations in this figure are the same as 
Fig. 1. 



C. &csl 

The OPE for this current is given as follows 

^i-nc)tf) = _^£,i„_4_^«(„-i)-„j[-L(„)l=. 
^n^tf) = !^!i(«±^£*^«(„-i)-™jl-L(„)13, 

+mc(l - «)(3 - 5^)?^ + 2Mm^|[-L(ii)]2 
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FIG. 5: Similar figure as Fig. 2 for the current @csi with Sq = (3.3 GeV)^. Here each label represents 
contribution from each term in Eq.(28). 



-imnw g2 = ^^^y^^yw^^ / ^^i^ ^.^ .^j^ 

TT 9 ■ Jq 

-Imn»(g2) = MM + S{q' - ml). (28) 

TT 54 

Note here again that the superscripts correspond to the diagrams shown in Fig. 4. The 
dimension-5 condensate involving is related to the quark-gluon condensate via (qD'^q) = 
{qgcr ■ Gq)/2. Similar relation holds for the corresponding strange-quark condensate. The 
correction to this relation is proportional to square of the quark mass which should be very 
small even for the strange quark. Fig. (5) represents the OPE as defined in Eq.(25) with 
the imaginary part in Eq.(28). We have only included a few terms in the OPE to show how 
each term contributes differently to the sum rule. As can be seen from the figure, the line 
denoted as "OPE", which is sum of the power corrections only, are much larger than the 
perturbative contribution. Moreover, the gluon condensate from diquarks is only a small 
fraction of the large higher order correction. This suggests that the OPE are not convergent 
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and it is very unlikely that the diquark correlation will remain an important mechanism in 
this configuration. 

D. @cs2 

The OPE for this current is given as follows 
ilmn(«)(g') 

TT 

-Imn(^)(g2) 

TT 

-Imn(^)(g2) 

TT 

-Imn('^)(g2) 

TT 

-Imn(^)(g2) 

TT 

-Imn(^)(g2) 

TT 

ilmn(5)(g2) 

TT 

-Imn('*)(g2) 

TT 

-Imn»(?2) 

TT 

Again note that the OPE diagram for each label is shown in Fig. 4. Fig. (6) represents 
the OPE as defined in Eq.(25) with the imaginary part in Eq.(29). Again, we have only 
included a few terms in the OPE to show a general trend of each contribution. For the 
negative parity case, the OPE has large contributions with alternating signs. The situation 
is better for the positive parity case, but again, the power corrections alternate in signs. 

From all the previous analysis on the OPE for the charmed pentaquark with and without 
strangeness, we find that the one without strangeness with diquark structure are most 
reliable, and are dominated by gluon condensate coming from diquark correlation. It is 
interesting to note that this result is consistent with the Skyrme model calculation which 
predicts a bound state of pentaquarks in the nonstrange sector [15]. In the following, we 
will perform a more detailed analysis with the stable structure well represented by the 
interpolating current Oc2- 
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FIG. 6: Similar figure as Fig. 2 for the current @cs2 with Sq = (3.3 GeV)^. Here each label represents 
each term in Eq.(29). 

V. QCD SUM RULES AND ANALYSIS 

A. The couplings to the DN continuum, Adtv 

As discussed before, it is important to subtract out the contribution from the DN con- 
tinuum. For that, one needs to know the couphng strength Xdn- Here we determine this 
for the currents without strange quarks, Xdn,c2- In the case of 9"*" (1540) [37], the soft-kaon 
theorem was used to convert the external Kaon state, corresponding to the D meson states 
in Eq.(13) and Eq.(17), to a commutation relation of the operator and the corresponding 
axial charge. The strength of the resulting five-quark operator with an external nucleon 
state was then obtained from a separate nucleon sum rule analysis with the same five-quark 
nucleon current. However, applying the soft D meson limit will obviously not work in the 
present case. 

Instead, we determine the coupling strength directly from the sum rule method. To do 
that, we eliminate contribution from the low-lying pole by introducing the additional weight 
ly(g^) = — rriQ in Eq.(5). We will take me = 3 GeV, and confirmed that changing it by 
±200 MeV will have less than 5 % effect on the Xdn value. This way of eliminating a certain 
pole is sometime used in QCD sum rules [45, 46]. Then, substituting the corresponding 
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FIG. 7: The \\dn,c2\'^ from the sum rule for z = g (sohd) hne and i = 1 (dashed line). The upper 
(lower) sohd or dashed lines in this case are for = (3.8 GeV)^ [^o = (3.7 GeV)^]. 



imaginary parts, we find, 

f^dg^ e-^V^-(g^-m|)llmnr(g^) ^. ^ ^ 
\^dn\ = —Q , {i = i-,Q) (30) 

where the i = 1, g in ImH represent the part proportional to 1 or ^ in the respective imaginary 
part, and Imll^^ is the spectral density in Eq.(15) or in Eq.(19) without the |A£»jvP- 

Figure (7) shows the plot of Eq.(30). The two dotted (solid) hues represent boundary 
curves with the least Borel mass dependence for the Xdn from the 1 (g) sum rules. A^jv 
should not only be independent of the Borel mass but also independent of the sum rule 
from which it is obtained. However, the results coming from either i = q oi i = 1 sum rule 
differ slightly. Inspecting the OPE, one finds that the contributions from higher dimensional 
operators are consecutively suppressed for the i — q sum rule, while that is not so for the 
i = 1 sum rule. Therefore, the value from the former sum rule should be more reliable. 
Nonetheless, to allow for all variations, we will choose the following range for the |Ai5jvP 
values, 

2 X 10-^GeV^° < \Xdn,c2\^ < 3 x IQ-^GeV^^. (31) 
Similar attempts to determine Xdn,ci give vastly different values from either i — q or i — 1 
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FIG. 8: The left figure shows the left-hand side of Eq.(32) using @c2 for positive parity case with 
|ADiv,c2p = 2 X 10-5 GeV^° (dashed line) and \Xdn,c2\'^ = 3 x 10"^ GeV^° (dot-dashed line). The 

■,c2| 

different threshold parameters. See Eq. (8) for A and B. 



solid line is when there is no DN continuum, \X£)n,c2\'^ = 0. The right figure is obtained with 



sum rules. This reflects the non-convergence of OPE from which one can not expect a 
consistent result. 



B. Parity 



Wc will now concentrate on the sum rule obtained from 0^2. Using the dispersion relation 
in Eq.(lO) and the spectral density in Eq.(20), one finds the following sum rule, 



I A 



±,c2|^e" 



/M2 



dqo e 



Pope(?o) - Pdn{%] 



(32) 



As can be seen from Fig. (8), the left hand side of Eq.(32) is positive for positive parity 
case. For \\dn,c2\'^ = (the sohd lines), we have chosen the continuum threshold to 
be 3.4 GeV and 3.3 GeV, which gives the most stable pentaquark mass as we will show in 
the next subsection. Similar method was used to obtain the continuum thresholds when 
\^DN,c2\^ 7^ 0. However, Fig. (9) shows that the corresponding sum rule is negative for the 
negative parity case, suggesting that there can not be any negative parity state. This result 
also confirms the non-convergence of the OPE for the negative parity case, from which a 
consistent result can not be obtained. This can also be expected from the constituent quark 
picture. The two diquarks in Oc2 current have opposite parities and, when they are combined 
with the antiquark, the configuration should be dominated by the positive-parity part in the 
nonrelativistic limit. 
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FIG. 9: The left figure shows the left-hand side of Eq.(32) using @c2 for negative parity case with 
|AcAr,c2p = 2 X 10-5 GeV^° (dashed line) and \\dn,c2\'^ = 3 x 10"^ GeV^° (dot-dashed line). The 
solid line is when there is no DN continuum. The right figure is obtained with different threshold 
parameters. 



C. Mass 

The sum rule for the ©c mass is obtained by taking the derivative of Eq.(32) with respect 
to The solid and dashed lines in Fig. (10) represent the mass for two different \dn 

values. The threshold parameters were obtained to give the most stable mass within the 
Borel window plotted. One notes that the inclusion of the coupling to the DN continuum 
states, the mass reduces to smaller values to below 3 GeV. The curve with Xdn,c2 = 2 x 10"^ 
GeV^'^ lies between the solid and dashed lines in Fig. (10). This suggests the possibihty 
that the heavy pentaquark might actually be bound; namely, lies below the DN threshold. 
This is consistent with the constituent quark model picture, where one expects the diquark 
correlation to be more dominant than that of the quark-antiquark correlation as the partic- 
ipating antiquark becomes heavy. However, if this was the case, its existence can only be 
measured through its weak decay. 



VI. SUMMARY 



Wc have performed the OPE and QCD sum rule analysis for heavy pentaquark with and 
without strangeness with two different current each. We find that the OPE is convergent 
only for the non-strange pentaquark with diquark structure. The OPE for this structure is 
dominated by gluon condensate coming from diquark, which non-pcrturbatively represents 
their strong correlation. We find that the heavy pentaquark without strangeness has positive 
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FIG. 10: The mass obtained by taking the square root of the inverse ratio between left hand side 
of Eq.(32) and its derivative with respect to using 6c2- 

parity as reported earlier [34] and that its mass lies below 3 GeV, when the DN irreducible 

contribution is explicitly including in the phenomenological side of the sum rule. The picture 
that we described here does not work so well in the light pentaquark O"*", as the OPE are 
highly divergent [47] as can be seen in the picture of the OPE in the original sum rule paper 
for the light pentaquark state[29]. 
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